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We derive a general relation between the fine structure splitting (FSS) and the exciton polar- 
ization angle of self-assembled quantum dots (QDs) under uniaxial stress. We show that the FSS 
lower bound under external stress can be predicted by the exciton polarization angle and FSS un- 
der zero stress. The critical stress can also be determined by monitoring the change in exciton 
polarization angle. We confirm the theory by performing atomistic pseudopotential calculations for 
the InAs/GaAs QDs. The work provides a deep insight into the dots asymmetry and their optical 
properties, and a useful guide in selecting QDs with smallest FSS which are crucial in entangled 
photon sources applications. 
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The recent successful demonstration pj of entangled 
photon emission from a biexciton cascade process [3| in 
a single quantum dot (QD) represents a significant ad- 
vance in solid state quantum information applications. 
The "on-demand" QD entangled photon source has fun- 
damental advantages over traditional sources using spon- 
taneous parametric down conversion process which are 
probabilistic 3]- Furthermore, QD emitters can be 
driven electrically rather than optically, and therefore 
have many advantages in device applications 01 ■ How- 
ever, it is still a big challenge in the generation of high- 
quality entangled photon pairs from QDs. 

The key issue here is to suppress the fine structure 
splitting (FSS) [Fig. [T] (a)] of the monoexciton, which 
arises from the underlying asymmetry of the QDs. There 
have been many attempts to reduce the FSS in QDs, in- 
cluding thermal annealing choosing proper dot ma- 
trix materials @ , and growing the dots in high symmetry 
directions However, these methods can only reduce 
the FSS to the level of about 10 ^eV, which might still 
be too large to produce high-quality entangled photon 
pairs. It has been shown that the FSS can be effectively 
tuned by an external magnetic field [l[ and electric field 

Q . Perhaps a more convenient way to tune the FSS 



in a QD is via uniaxial stresses [10|. Singh et al [11 1 
showed that the FSS can be tuned to zero under uniaxial 
stress for an ideal QDs with C2v symmetry. However, for 
a general QD, which has Ci symmetry, there is a lower 
bound for the FSS when an external stress is applied. It 
was unclear which or what kind of QDs may have the 
smallest FSS lower bound under stress. Therefore, it is 
still an open question in selecting QDs that are suitable 
for entangled photon emitters. 

The purpose of the letter is to establish a general rela- 
tionship between the asymmetry in QDs, the exciton po- 
larization angle and the FSS under uniaxial stress, there- 
fore provides a useful guide in selecting QDs that have the 
smallest FSS lower bounds for applications such as entan- 
gled photon emitters. We show that QDs in which the 
exciton has polarization closely aligned along the [110] 
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FIG. 1; (Color onfine) (a) Four lowest energy levels of mo- 
noexciton in a QD, including two dark states (broken lines) 
and two bright states (solid lines). The energy difference be- 
tween the two bright states defines the FSS. (b) Polarization 
of exciton emissions under different symmetry. The black line 
corresponds to H2v without spin-orbit coupling, the blue line 
corresponds to H2v with spin-orbital coupling, and the red 
dotted line corresponds to C'l symmetry. 



or [100] directions have the smallest FSS under stress. 
The critical stress can also be determined by monitoring 
the change in the exciton polarization angle. The the- 
ory is further confirmed through an atomistic empirical 
pseudopotential calculations for InAs/GaAs QDs. 

We start from a general understanding of the relation 
between the QD symmetry and exciton polarization an- 
gles shown in Fig. [ijb). An ideal InAs/GaAs QD has C2v 
symmetry. Without spin-orbit coupling, the two bright 
exciton states belong to two different irreducible repre- 
sentations |r2) and |r4). The polarizations of the two 
bright states are along the [100] and [010] directions re- 
spectively. When the spin-orbital interaction is included, 
the two bright states belong to irreducible representa- 
tion ]r2) ± ilr4) [l3], therefore the polarizations of the 
emission lines should be exactly along the [110] and [110] 
directions. For a general dot, the symmetry is further 
lowered to Ci owing to structural imperfections or alloy 
randomness [isl ]. the polarization angle will depart from 
the [110] and [110] directions, i.e., 9 ^0. 
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When uniaxial stress is applied, the exciton Hamilto- 
nian can be written as, 



(a) p = o 



(b) a = 



H{n,p)=H2v + Vi+Vs{n)p, 



(1) 



where n is the external stress direction, and p is the 
magnitude of the stress. represent the Hamilto- 

nian of an idea QDs with C2v symmetry, whereas Vi 
lower the dot symmetry to Ci, due to local structure 
deformations, alloy distribution 13| and interfacial ef- 
fects 14] etc. Vs{n)p is the potential change due to 
the external stress. We neglect the higher-order 0{p^) 
terms. This is justified by atomistic pseudopotential cal- 
culations, which show that these terms are negligible up 
to 200 MPa. The eigenvectors of the two bright states 
of iJ2, are |3) = {T^} + t\r^) and |4) = jPs) - i\Ti), 
with corresponding eigenvalues £"3 and E4, respectively. 
The energy levels are schematically shown in Fig. [TJa) 
in sohd lines. The difference A — lE^ — E^] is the FSS. 
The other two states |1) and |2) are optically dark, also 
shown in Fig. [IJa) in broken lines. In the absence of an 
in-plane magnetic field, the coupling between the dark 
states and bright states is negligible. We therefore write 
the Hamiltonian in the space spanned by the two bright 
states. 



H 



E + 5 + azp _ K + Pp 
K + /3p E — S + a^p 



(2) 



where ^+,5 = (3|i/2.+l^i|3), E^5 = (4|i/2.+l^i|4). a, = 
(i|14(n)|i) (i=3, 4), n = (3|1/i|4) and /? = (3|K(n)|4). 
Because the Hamiltonian has a time-reversal symmetry, 
all parameters can therefore be set to real values for sim- 
plicity. Diagonalization of the Hamiltonian yields eigen- 
values, 



- , P7 ± \/4(/3p + + {ap + 25Y 
E±=E-\ , (3) 

where a — 03 — 014 and 7 = 03-1-04. 
The eigenvectors of the two states are \tp±) = 

- "''^"^^.X+rr^' ^l^) + |4>- Since aU param- 
eters are real, the two states are linearly polarized [l5| . 
We calculate A(p) for QDs under uniaxial stress p, 



A{p) = y^AiPp + k)2 + {ap + 2<5)2 



(4) 



Because is a quadratic function of p, the lower bound 
of FSS and the corresponding critical stress can be cal- 
culated analytically. 
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The polarization angle vs p can be calculated from. 
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FIG. 2: (Color online) FSS vs p for (a) /3 = and (b) a = 0. 
The circles above show the exciton polarization angles at pc 
and far away from pc- 



At the critical stress Pc, 
where a is the 
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Obviously, 9 changes with p 
we have tan (6'^) 
sign of (—2/3(5-1- an). It is interesting to find that the 
polarization angle at the critical stress is independent of 
5 and k, but only on the ratio of /3/o. Thus the values 
of o, /3, K, and 5 can be uniquely determined using the 
relationship between A and p, and the polarization angle 
at pt,. 

For a dot with C2V symmetry, the external stress along 
the [110] and [110] directions would not change the sym- 
metry of the dot. For a general dot, the symmetry of 
the dot is lowered to Ci due to alloy distribution and 
structural asymmetry [ij, [3], which change little un- 
der external stress. Therefore, for stress along the [110] 
and [110] directions, we have /3 wO. Another special case 
arises when stress is applied along the [100] direction, for 
which we have o=0 by symmetry. The results for these 
two special cases are schematically shown in Fig. [2] (a), 
(b) respectively. In Fig. [H^a), the change in FSS with p 
is determined by ]o], and the lower bound is determined 
by 2]k] at Pc — ~2d/a. At p = pc, the polarizations are 
along the [100] or [010] directions whereas far away from 
Pc, the polarizations are along the [110] or [110] direc- 
tions. The results for o = shown in Fig. [5] (b) are 
totally different. In this case, the change in FSS with 
p is determined by 2]/3], and the lower bound is 2\5\ at 
Pc = ~n/ p. At p — Pc, the polarization is along the 
[110] and [110] directions, whereas far away from pc, the 
polarization is rotated into the [100] and [010] directions. 
Therefore, the polarization angle of the emission lines can 
be used to determine the critical point pc in experiments. 
The above picture is also correct for QDs with C2V sym- 
metry, where n — and the results in Fi g. [2] (a) are then 
reduced to the results presented in Ref. [111 ]. 

The theory also provides a simple way to determine the 
FSS lower bound of a QD before applying the external 
At p = 0, we have Ao = 2\/P 
s 



(6) 



stress. At p = 
tion angle tan(6') 
that K = — Aq 



sin(26')/2 and 6 = Aq 



-f k2 and polariza- 
It is easy to show 
cos(26l)/2. Thus 
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TABLE I: Parameters for Ina;Gai_a;As/GaAs QDs under uni- 
axial stress, n is the direction of external stress, a and /? are 
in units of ^eV/MPa, k and S are in the units of /ieV, the 
critical stress pc is in MPa and the lower bounds Sb is in units 
of /xeV. The unit for base diameter D and height h of the dots 
is nm. 



QDs 




a 


/3 


5 




Pc 


5b 


Lens (a::=0) 


[110] 


0.36 





-8.46 





47 





73=20, h=3.5 


[110] 


-0.37 





-8.55 





-46 







[100] 





-0.095 


-8.65 








17.3 


Lens (x = 0.6) 


[110] 


0.14 





-3.36 


-2.90 


49 


5.8 


73=25, /i=3.5 


[110] 


-0.14 





-3.36 


-2.86 


-47 


5.7 




[100] 





-0.047 


-3.38 


-2.90 


-62 


6.8 


Pyramid (a;=0.6) 


[110] 


0.13 





-1.32 


-0.64 


21 


1.3 


7)=25, /i=3.5 


[110] 


-0.13 





-1.41 


-0.61 


-22 


1.2 




[100] 





-0.048 


-1.32 


-0.69 


14 


2.6 


Elongated (2;=0.6) [110] 


0.14 





2.76 


-0.72 


-40 


1.4 


I^[iioi=26, 


[110] 


-0.14 





2.73 


-0.78 


38 


1.6 


7)[no]=20, /i=3.5 


[100] 





-0.051 


2.62 


-0.92 


-18 
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by measuring the polarization angle and FSS at p=0, we 
can uniquely determine the values of 2|k| and 2\6\ which 
are the FSS lower bounds for the stress along the [110] 
([110]) and [100] directions. For entangled photon source 
applications, \6\ or |k| need to be smaller than 1 fieY. 
Therefore our results provide a useful guide in selecting 
QDs for entangled photon sources, i.e., one should pick 
QDs for which the polarization angle is as closely aligned 
as possible to the [110] or [100] direction. 

To confirm the above analysis, we perform numeri- 
cal calculations of the FSS and polarization angle of 
InAs/GaAs QDs using an empirical pseudopotential 
method 0, Il3|- We model the InAs/GaAs QDs by em- 
bedding the InAs dots into a 60x60x60 8- atom GaAs su- 
percell. To calculate the exciton energies and their FSS, 
we first have to obtain the single-particle energy levels 
and wavefunctions by solving the Schrodinger equation. 



V'i(r) = a ^pi{r) 



(7) 



where Vps{r) 



so - 



'Y.^^aVa{r-'Ri^a)+i^piczo IS the to- 
tal pseudopotential and ^^(r — Rj q) is the local screened 
atomic potential at the equilibrium atom position R^^a 
obtained by minimizing the total strain energies under 
the given stress 11| using the valence force field method 
[l8| . Vso is the spin-orbit interaction, and Wpiozo is the 
piezoelectric potential The single particle energy 

levels are calculated using a linear combination of bulk 
bands method [lB|- The exciton energies are then cal- 
culated via many-particle configuration interaction (CI) 
method [20| . in which the exciton wavefunctions are ex- 
panded in Slater determinants constructed from all con- 
fined electron and hole single-particle states. 

The change in FSS under external stress is purely an 
atomistic effect, because the macroscopic shapes of the 
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FIG. 3: (Color online) Left panel: The FSS as a function of p 
for (a) pure lens-shaped, (c) alloy lens-shaped, (e) alloy pyra- 
midal and (g) alloy elongated InAs/GaAs QDs. Right panel: 
Corresponding exciton polarization angle as a function of p. 
The open squares (circles) are the results for stress along the 
[110] ([100]) direction calculated from the atomistic pseudopo- 
tential method, whereas the solid lines are the results of the 
theory. Solid squares and circles indicate critical stresses. 



dots change little (less than 0.1%) under such stresses and 
should not affect the FSS. Furthermore, we find that in- 
cluding piezoelectricity or not gives essentially the same 
results, suggesting that piezoelectricity is not responsi- 
ble for the FSS change. In fact, the change in FSS is 
due to changes in the underlying atomic structure. We 
have calculated the FSS of more than 13 dots under ex- 
ternal stresses along the [110], [lIO] and [100] directions. 
The behaviors of FSS are almost symmetric for stresses 
along the [110] and [110] directions, i.e., the effects of 
tensile stress along the [110] direction is almost identical 
to the effects of compression along the [110] direction. 
The results for some typical dots are shown in Fig. [3l 
whereas the geometry and other parameters of these dots 
are listed in Table. HI In the left panels of Fig. [3l we plot 
the FSS vs p along the [110] (black square) and [100] (red 
circle) directions. The solid lines are fitted from theory 
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using Eq. (|4]). The right panels show the correspond- 
ing polarization angle 9 vs p where the solid lines are the 
theoretical predictions using Eq. We use /3=0 (a^O) 
for p along [110] ([100]) direction. As can be seen, the 
agreement between numerical calculations and theory is 
remarkable. 

(1) Pure InAs/GaAs QDs with C2v symmetry. In Fig. 
[3] (a) we show the FSS vs p along different directions for 
a pure lens-shaped QDs with base D =20 nm and height 
/i=3.5 nm. When the stresses are directed along the [110] 
and [110] directions, the FSS can be tuned exactly to zero 
[Tlj . The polarization angle 9 is constant (90°) below pc 
— 47 MPa and jump to 180° after pc as shown in Fig. 
[3](b). However, if the stress is along the [100] direction, 
the FSS can not be tuned to zero, in agreement with 
previous results ll|. The polarization angle rotates tol- 
as seen in Fig. [3] (b). At p — Pc, the 
We also calculate pure pyrami- 



lowing Eq. 

polarization angle is 90° 
dal and elongated QDs, which do not have macroscopic 
cylindrical symmetry but still retain C2v symmetry, and 
find similar features. 

(2) Alloy InQ ^GaoAAs/GaAs QDs with Ci symmetry. 
For alloy dots, the symmetry is lowered to Ci . The FSS 
has lower bound under the uniaxial stress [ll|, as shown 
in Fig. [3](c) (e) and (g), for different dot geometries and 
sizes. The corresponding parameters are summarized in 
Table m The stress dependence of the polarization angles 
is also in excellent agreement with theory as shown in Fig. 
[3](d), (f) and (h) for the three dots. All three dots have 
polarization angle 6'c=135° at Pc if the stress is along the 
[110] direction, and ^£=90° (or 0°) if the stress is along 
the [100] direction, as predicted by the theory. Among 
the three alloy dots, the lens-shaped QD [Fig. [21(c)] has 
the largest lower bounds ^ 7 /zeV at pc ^ 62 MPa along 
the [100] direction. At p—Q, the lens-shaped QD has 9 
= 110°, compared to 6* = 103° for the pyramidal dot 
and 169° for the elongated dot. The polarization angle 
of the lens-shaped QD deviates from the [110] (or [110]) 
direction most, and hence has the largest FSS at Pc, as 
predicted by the theory. We obtain similar results when 
the stress in along the [100] direction. 

In the calculations, we find that a is not very sensitive 
to the QD shape, but changes with alloy compositions. 
For example, for pure dots, \a\ ^ 0.2 - 0.4 /xeV/MPa for 
p along the [110] direction, whereas for In composition 
x—0.6, \a\ reduces to 0.1 - 0.2 /ieV/MPa. (3 also has sim- 
ilar features (for p along the [100] direction), with \(3\ ~ 
0.05 - 0.1 /ieV/MPa for pure dots, and \l3\ - 0.04 - 0.05 
/ieV/MPa for alloy dots with a::=0.6. We also calculate 
alloy dots of the same geometry but with different alloy 
distributions [lH and find that the alloy distribution does 
not significantly change the values of a and /?. In con- 
trast, 6, K, and the polarization angle 9 at p—0, change 
dramatically from dot to dot, in agreement with recent 
experiments fSl']. However, in all cases, the behaviors of 
the FSS and polarization angle under stress are in excel- 



lent agreement with our theoretical predictions. 

To conclude, we have established a general relation- 
ship between the asymmetry in QDs, the exciton polar- 
ization angle, and the FSS under uniaxial stress. We 
showed that the FSS lower bound under external stress 
can be predicted by the polarization angle and FSS under 
zero stress. The critical stress can also be determined by 
monitoring the change in exciton polarization angle. The 
work therefore provides a useful guide in selecting QDs 
with smallest FSS which is crucial for entangled photon 
sources applications. 

The authors thank A. J. Bennett for bringing Ref. 
to our attention. LH acknowledges the support from 
the Chinese National Fundamental Research Program 
2011CB921200 and National Natural Science Funds for 
Distinguished Young Scholars. 

Note added: After submitted the paper, we became 
aware of Ref. i. There, the exciton FSS of InAs/GaAs 
QDs is tuned via an electric field along the [001] direction, 
which has the same symmetry as applying stress along 
the [110] and [1-10] directions (i.e., the case of P=0). 
Determining the degree of agreement between the present 
theory and the experiment of Ref. is a promising avenue 
for future research. 
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